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We study the frequency and angular dependences of Cherenkov radio pulses originated by the 
excess of electrons in electromagnetic showers in different dense media. We develop a simple model 
to relate the main characteristics of the electric field spectrum to the properties of the shower such 
as longitudinal and lateral development. This model allows us to establish the scaling of the electric 
field spectrum with the properties of the medium such as density, radiation length, Moliere radius, 
critical energy and refraction index. We normalize the predictions of the scaling relations to the 
numerical results obtained in our own developed GEANT4-based Monte Carlo simulation, and we 
give a unified parameterization of the frequency spectrum and angular distribution of the electric 
field in ice, salt, and the lunar regolith, in terms of the relevant properties of the media. Our 
parameterizations are valid for electromagnetic showers below the energy at which the Landau- 
Pomeranchuk-Migdal effect starts to be relevant in these media. They also provide an approximate 
estimate of radio emission in hadronic showers induced by high energy cosmic rays or neutrinos. 
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I. INTRODUCTION 

Neutrinos are the most suitable candidates to extend 
astronomical observations to the ultra high energy (UHE) 
regime, above ~ 1 TeV. Unlike charged particles, they 
point directly to the source where they were produced. 
Unlike neutrons, they are stable. Unlike photons, they 
can penetrate large amounts of matter, and they are not 
attenuated by cosmic backgrounds. They are not optimal 
candidates though, due to their small interaction cross 
section, and to their small expected fluxes at ultra high 
energy However these two latter difficulties may 

be overcome with a large enough detector volume of at 
least 1 km^ Q. 

The “conventional” technique to achieve a large detec¬ 
tor volume exploits the observation of Cherenkov light 
emitted by high energy neutrino-induced muons and 
showers in water and ice iiHii An alternative to 
it is the search for coherent Cherenkov radio pulses in the 
MHz-GHz radio frequency range, produced by neutrino 
induced showers in dielectric, transparent, dense media. 
When the wavelength of the emitted radiation is larger 
than the typical dimensions of the shower the emission is 
coherent. The contribution to the electric field from pos¬ 
itive and negative particles would approximately cancel 
out were it not for the existence of an excess of electrons 
over positrons in the shower M- The charge asymmetry 
arises from photon and electron scattering on the sur¬ 
rounding medium that sweeps electrons into the shower, 
and from positron annihilation in flight that contributes 
to the excess charge by terminating positron trajecto¬ 
ries. Due to the coherent nature of the radio emission 
the power in radio waves scales as the square of shower 
energy. 

Dense media have the advantage that shower dimen¬ 
sions are of the order of meters, and coherence extends 
to high frequencies, typically ~ GHz, where more power 


is available because the coherently radiated electric field 
is proportional to frequency. Besides, large formations 
of dense, transparent media such as ice and salt exist 
in nature having attenuation l eng ths in radio frequen¬ 
cies of a few hundred meters jlll IT^ IT^ . They have 
the advantage of being located in low noise environments 
such as the South Pole or deep salt mines. A single an¬ 
tenna in such a medium makes a detector of effective 
volume of ^ 0.1 km^ water equivalent HHH The cost- 
effectiveness of antennas also adds to the attractiveness 
of the radio technique. 

Several experiments are searching for neutrinos ex¬ 
ploiting the radio technique in dense media. The RICE 
experiment M is an array of antennas buried in the 
Antarctic ice cap that has imposed stringent limits on 
astrophysical neutrino fluxes in the energy range above 
~ eV Ea The Goldstone Lunar Experiment 

(GLUE) has used two large radiotelescopes to search 
for pulses produced by neutrino-induced showers on the 
Moon’s regolith 0 . With no neutrino candidates in 
120 hours of observation, GLUE has established compet¬ 
itive upper bounds on astrophysical neutrinos of energy 
above ^ 10^® eV 0. The technique could be of inter¬ 
est to search for ultrahigh energy cosmic rays 00 
and is being explored using the Westerbork array of ra¬ 
diotelescopes in the Netherlands 0 |. It will also be ex¬ 
ploited^^ the proposed LUNASKA experiment in Aus¬ 
tralia 01. A recently approved experiment is ANITA 
1^. a balloon-borne antenna that circles the Antarctic 
continent scanning for large radio pulses. The idea of us¬ 
ing large salt domes as targets for neutrino interactions 
has also been explored and there is a proposed experi¬ 
ment named SalSA [0, as well as other initiatives such 
as the ZESANA proposal in the Netherlands 0 . 

It is also very important to remark that the dominant 
emission mechanism in dense media has been experimen¬ 
tally confirmed in accelerator measurements at SLAG, 
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using bunches of bremsstrahlung photons as projectiles, 
and sand and salt as target medium 1^ . This ex¬ 
periment has tested the theoretica.1 predictions originally 
made by Askar’yan in the 1960’s flfl . 

Given the current and future experimental efforts, the 
motivation and aim of our work are clear: A good un¬ 
derstanding and comprehensive characterization of the 
dependence of the radio signal on frequency and obser¬ 
vation angle in different dense media are needed to inter¬ 
pret the data that is being collected mui , as well as 
to evaluate the capabilities and potential of the planned 
experiments looking for neutrinos and cosmic rays. 

In this paper we show that many of the observables of 
the electric field spectrum emitted by an electromagnetic 
shower developing in dense media, scale to a few percent 
level with several properties of the media, such as density 
(p), radiation length (Xq), Moliere radius (Rm), critical 
energy (Ec), and index of refraction (n). This scaling 
allows us to predict the properties of radio pulses in dif¬ 
ferent media without performing time-consuming Monte 
Carlo simulations. Our work provides further insight into 
the close relation between shower development and radio 
emission, which may help designing future experiments 
exploiting the radio technique, as well as assessing the 
potential of those currently under construction. 

This paper is structured as follows. In Section UTI we 
briefly review radio emission from electromagnetic show¬ 
ers. In Section fllll bv means of a simple toy model we re¬ 
late the main features of the radiopulse to general charac¬ 
teristics of shower development. This allows us to obtain 
the scaling of radiopulse with the relevant parameters 
of the medium. In Section im we normalize these scal¬ 
ing relations using our own developed GEANT4-based 
Monte Carlo simulation to numerically calculate 

the electric field spectrum in ice, salt and the lunar re- 
golith. We also give in this section unified parameteri- 
zations of our numerical results to be used in practical 
applications and establish their range of applicability. 


II. THE RADIO PULSE SPECTRUM 

We provide here a short review on Cherenkov radio 
emission from electromagnetic showers for the sake of 
completitude. Further information can be found in ref¬ 
erences 

The Fourier time-transform in the Fraunhofer limit 
of the electric field radiated by a charge q moving at 
constant velocity v for an infinitesimal time interval 
-I- 6t) starting at the position fq is given by [^: 


Eiu;,r) = ( 1 ) 

27reoC^ r y 

with y = TTvStll — nP cos 0). Here eo (c) is the permittiv¬ 
ity (speed of light) in vacuum, (n) is the relative mag¬ 
netic permeability (the refraction index) of the medium, 


Lo the angular frequency, fc is a unit vector in the direc¬ 
tion of observation, and v± is the projection of v in the 
direction perpendicular to k. r is the distance from the 
charge to the observer’s position. It is straightforward to 
deduce Fq. m from the inhomogeneous Maxwell’s equa¬ 
tions in the transverse gauge for a point charge current 
keeping only the R~^ radiation term. 

Fq. Q displays the most important characteristics of 
the radiation emitted by a charged particle travelling 
along a straight trajectory at a speed larger than the 
speed of light in the medium: 


1. The electric field spectrum rises linearly with fre¬ 
quency, and has a siny/y term, where y oc vSpl — 
n/3 cos 9) that corresponds to an angular diffraction 
pattern peaked at the Cherenkov angle 9c given by 
cos 0c = (n/3)“^. The first zeros of the diffraction 
pattern happen at. 


vSt^vPp'^ — I ’ 

that is the width of the peak around the Cherenkov 
angle is proportional to the inverse of the length 
of the particle’s track, and to the inverse of the 
frequency of observation, as corresponds to the 
diffraction pattern produced by a single slit. For 
observation at the Cherenkov angle there is no 
phase factor associated to the position of the parti¬ 
cle along the track and the maximum in the diffrac¬ 
tion pattern is achieved |4^. 

2. The electric field is proportional to the length of 
the particle track perpendicular to the direction of 
observation, i.e. it scales as w<5<sin0. 

3. The field is 100% polarized in the plane containing 
the observation direction and the particle trajec¬ 
tory. 

In order to calculate the electric field produced by a 
shower, the contributions from each particle track, as 
given by Fq. CJ, have to be summed taking into ac¬ 
count their relative phases which stem from their dif¬ 
ferent spacetime positions and different /3. To do this 
numerically, the tracks can be subdivided in small steps 
so that V is approximately constant in each step. 

The contribution to the total electric held from the 
particle tracks in a shower is strongly dependent on ob¬ 
servation angle and frequency. The problem resembles 
a classical problem of interference in which the source 
emits in a coherent fashion as long as the wavelength of 
observation is insensitive to (i.e. larger than) the hne de¬ 
tails of the emitting region. Further information on the 
characteristics of the electric held can be found in refer¬ 
ences [^13 EE US HE 113 • In Fig. n]we show the results 
of our own developed GFANT4-based simulation of the 
electric held spectrum emitted by a 10 TeV shower in ice 
|3(T|. The spectrum is plotted for different observation 
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FIG. 1: Frequency spectrum of the electric field in ice ob¬ 
served at several angles from the shower axis: 8 = 8c — 55.8° 
(crosses); 8c + 10° (bold squares); 8c — 10° (empty squares); 
0c+ 20° (bold circles); 8c —20° (empty circles); 0c+ 30° (bold 
triangles); 0c — 30° (empty triangles). The results were ob¬ 
tained using our own developed GEANT4-based Monte Garlo 
simulation. The lines joining the points are just to guide the 
eye. 


angles 0 with respect to shower axis. The electric field 
spectrum rises linearly with frequency up to a turnover 
point which depends on the observation angle. 


III. TOY MODEL OF SHOWER 
DEVELOPMENT 

In this section we describe a simple toy model of shower 
development that allows us to relate the main features of 
the electric field spectrum to the properties of the shower 
producing the field. 

When a photon or electron of energy E incides on 
a thick absorber it initiates an electromagnetic shower. 
The number of electrons, positrons and photons in the 
shower rapidly grows due to bremsstrahlung and pair 
production which are the dominant processes at high 
energy. The energy of the primary particle initiating 
the shower quickly degrades because it is shared among 
an increasing number of secondaries, until their energy 
reach the so-called critical energy i?c, at which other 
processes become important, and electrons are as likely 
to lose all their energy through ionization as to radi¬ 
ate hard bremsstrahlung photons. As a result of the 
shower development, the number of particles increases 
up to a maximum fVmax, and exponentially decreases af¬ 
ter it. Using the Heitler model ^3 that assumes that 
in each interaction the energy is shared equally between 
the two secondaries, it is straightforward to show that 


.shower axis 
i 1 



>'.' 

FIG. 2: Toy model of a shower. L is the effective length of 
the shower proportional to the radiation length Lq and R is 
the shower width, proportional to the Moliere radius Rm- 

Nmax = 2P = E/Ec, where p is the number of interac¬ 
tions before the shower reaches its maximum. Further¬ 
more, in this model electrons and photons are assumed 
to travel a radiation length Lq = X^jp before interact¬ 
ing. As a consequence the total tracklength T due to 
the charged particles in a photon initiated shower can be 
approximated by | 23 , 

T-(2 + 22 + ... + 2P)Lo-A'maxTo = ^ —. (3) 

Ec p 

The particles in the shower also spread in the trans¬ 
verse dimension - perpendicular to shower axis - mainly 
due to multiple scattering. We can view a shower as a 
thin pancake of particles travelling at the speed of light 
{(3 = 1) along a length L = kLLo = k^Xo/p proportional 
to the radiation length (Fig.|21). The pancake has a lateral 
width R = kuRMlp proportional to the Moliere radius 
Rm , the typical scale of the lateral spread of the shower 
(Fig. 13 . Both k]^ and k^ are normalization constants 
of the model to be determined in numerical simulations 
(see below). 

This very simple model of shower development, as 
shown below, allows us to establish the relation between 
the frequency, angular distribution and absolute mag¬ 
nitude of the electric field, and the parameters of the 
medium in which it is emitted. 


A. Electric field spectrum: shape 

The overall time duration of the radio pulse for a given 
observation angle from shower axis (0) is an exceedingly 
important parameter. It is related to the turnover fre¬ 
quency at which the spectrum departs from a linear be¬ 
havior with frequency. There are two relevant time scales 
in our simple toy model, namely, dfn the time delay be¬ 
tween two light rays emitted at the two end points in the 
longitudinal development of the shower (points 1 and 2 
in Fig. 12), and dlR, the time delay between two rays emit¬ 
ted from the two ends of the lateral spread of the shower 
(points 2 and 3 in Fig. [2). The dominance of one or the 
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other delay - i.e. the time duration of the pulse - de¬ 
pends on the observation angle. It is straightforward to 
obtain that: 

= — (1 ~ ncos0) = — (cos0c — COS0), (4) 

c c p 

where we have used cos 0c = ^/n, and 

otR = —— sin0 =-sin0. (5) 

c/n c p 


Here is a normalization constant to be determined 
in Monte Carlo simulations. Its numerical value is ex¬ 
pected to be different from despite the fact that both 
the turnover frequency of the spectrum away from the 
Cherenkov angle, and the angular distribution of the 
pulse are mainly determined by the longitudinal devel¬ 
opment of the shower. 

B. Electric field spectrum: magnitude 


For observation at the Cherenkov angle all the rays emit¬ 
ted along L at the same radial distance to the shower axis 
arrive at the same time at the observer, i.e. dtc vanishes 
at 0 = 0c. Only the time delay associated to the lateral 
spread of the shower remains and the electric field spec¬ 
trum increases linearly with frequency up to a turnover 
frequency given by 


i^r( 0 = 0c) ^ 


P 


1 


StR kuRu — 1 


( 6 ) 


where we have used sin 0c = V1 — cos^ 0c = Vn'^ — 1/n. 

The corresponding turnover frequency associated to 

(5tL is, 


t'L(0) 


1 


P _ c 

fciXo |1 —ncos0| 


(7) 


For observation away from the Cherenkov angle both 
6t^{d) and 6ti^{d) do not vanish. Eq. 0) illustrates the 
fact that 6th (0), the effective time duration of the longi¬ 
tudinal spread of the shower as viewed by the observer, 
grows from zero as the observation angle departs from 0c, 
and equals StR for observation angles 0+ = 0c -I- 69+ and 
0_ = 9c — 69- given by the condition 6th{9±) = 6thi{9±). 
For small values of 69±, 


69± 


_ knRM sin 9c _ 

khXo sin 0c T kuR+i cos 0c' 


( 8 ) 


As a result, if 0 G (0_, 0+), the shower electric field loses 
its full coherence at a turnover frequency z/r determined 
mainly by the lateral spread of the shower and given by 
an equation similar to Eq. © but for an arbitrary 0. 
Correspondingly, if 0 falls outside the (0_,0+) interval, 
the turnover frequency in the spectrum is determined 
by the longitudinal shower development, and is given by 
Eq. 0 . As a limiting case of Eq. 0 , in a shower having 
L >> R, 69± ~ R/L « 1 and the interference from 
different stages in the longitudinal shower development 
dominates the spectrum at essentially all angles, as ex¬ 
pected in a one dimensional shower where L is the only 
relevant scale. 

The width of the angular distribution of the electric 
field around the Cherenkov angle is given by an equation 
similar to Eq. © with an effective 6t which is also pro¬ 
portional to L. As a result it is easy to infer that it scales 
as 


A0 = 


P 


I 


/caAIq ^/n'^ — 1 


(9) 


The absolute normalization of the field at low frequen¬ 
cies - in the fully coherent region where the electric field 
increases linearly with frequency - is known to be de¬ 
termined by the distance travelled by the excess charge 
along the perpendicular to the direction of observation. 
Using the expression for the tracklength in Eq. ©, the 
absolute normalization of the electric field for an obser¬ 
vation angle 0 is expected to scale as 

-♦ E X 

r\E\ ^ kE V T sin0 ~ kE — —- v sin0, (10) 

Ec p 

where is a normalization constant to be determined 
in numerical simulations (see below), and r is the dis¬ 
tance from the source to the observer. It is important 
to notice that we are implicitely assuming that particles 
travel parallel to shower axis. As a consequence, the pro¬ 
jection of the tracks in the perpendicular to the direction 
of observation is simply accounted for by a factor sin 0 in 
Eq. m- 

IV. RADIO PULSE SPECTRUM IN 
DIFFERENT MEDIA 

A. Radio pulse simulations in dense media 

Using our own developed GEANT4-based Monte Carlo 
code 1^ , we have simulated electromagnetic show¬ 
ers in ice, salt and the lunar regolith. These media are 
suitable targets in which neutrino interactions can be ef¬ 
ficiently detected. It is worth remarking that we have 
checked these simulations against the ZHS Monte Carlo 

S , an independent code whose results were shown in 
to agree with those of CEANT4 to a few percent 
level in ice. Table Q] lists the relevant properties of the 
three media adopted in this work. 

In Fig. 0 we show the frequency spectrum of the elec¬ 
tric field emitted by a 10 TeV electromagnetic shower 
observed at the Cherenkov angle and at 0 = 90 deg. 
in ice, salt and the Moon’s regolith, as obtained in our 
CEANT4-based simulations. In Fig. 01 we plot the an¬ 
gular distribution of the electric field at 100 MHz and 1 
CHz. It is obvious by inspection that the shape of the fre¬ 
quency spectrum is fairly universal, only the normaliza¬ 
tion and the turnover frequency in the spectrum change 
from medium to medium and scale with its properties, 
as will be quantitatively shown in the next subsection. 
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TABLE I: Some relevant properties of the materials consid¬ 
ered in this work: density p, radiation length Xq, Moliere 
radins Rm, critical energy Ec, and index of refraction at ra¬ 
dio frequencies n. Last column is the energy above which the 
Landau-Pomeranchuk-Migdal effect starts to affect shower de¬ 
velopment using the definition in m. 


Medium 

P 

[g cm“®] 

Xo 

[g cm-2] 

Rm 

[g cm-2] 

Ec 

[MeV] 

n 

FIlpm 

[TeV] 

Ice 

0.92 

36.08 

10.35 

73.0 

1.78 

~ 2400 

Salt 

2.05 

22.16 

12.09 

38.5 

2.45 

~ 660 

Moon 

1.80 

22.59 

11.70 

40.0 

1.73 

~ 770 



FIG. 3: The frequency spectrum of the electric field emitted 
by a 10 TeV shower in the direction of the Cherenkov angle 
and at 90 degrees in different dense media. Empty squares: 
ice, 9 = 9c', empty circles: salt, 9 = 9c; empty triangles: lunar 
regolith, 9 = 9c- Filled squares: ice, 0 = 90 deg. ; filled 
circles: salt, 9 = 90 deg. ; filled triangles: lunar regolith, 
9 = 90 deg. The lines are the fits of the Monte Carlo results 
to Eqs. dJ and d, using the values of the parameters in 

Table El 


B. A unified parameterization 

From Eqs. ®, 0 , 0 and m, we have obtained 
a unified parameterization of the electric field spectrum 
and its angular distribution around the Cherenkov an¬ 
gle in terms of the relevant properties of the media that 
control radio emission. The electric field spectrum at the 
Cherenkov angle can be parameterized by the expression: 


r\Eii')\ec 


E Xo „ 1 

-Sint^c -;—^ r 

Ec p 1 MHz 1 + {vIvr) 


. . 

where vr is given in Eq. ®. The parameterization is 
a good approximation to the electric field spectrum for 
frequencies typically below u =10 GHz where coherent 
effects are still very important. 

Away from the Cherenkov angle (i.e. for 0 outside 



FIG. 4: Angular distribution of the electric field emitted by 
a 10 TeV shower in ice, salt and the lunar regolith. The 
distributions are shown at i/ = 100 MHz and v = 1 GHz. 
Empty squares: ice, zz = 1 GHz; filled squares: ice, v = 100 
MHz. Empty circles: salt, v = 1 GHz; filled circles: salt, 
V = 100 MHz. Empty triangles: lunar regolith, u = 1 GHz; 
filled triangles: lunar regolith, v = 100 MHz. The lines are 
the fits of the Monte Carlo results to Eq. d using the values 
of the parameters in Table E] 


the interval {9-,0+) ) the electric field spectrum is well 
described by: 


r\E{v)\g 


Ec p I MHz 1 + {v/vLf 


( 12 ) 


where vl is given in Eq. 0 . This expression is generally 
valid up to frequencies ^ 2vr. For larger frequencies, 
the fine structure of the shower is relevant (even at the 
individual particle level), and clearly our simplified toy 
model does not account for it. 

The angular distribution of the pulse around the 
Cherenkov angle can be parameterized by a Gaussian 
peak modulated by a sin0 which accounts for the angu¬ 
lar behavior in the fully coherent region: 


T\E{e)\=r\E{v)\e^ 


sin 6* 
sin 9 c 


exp 


e-9c\ 

A9 ) ’ 


(13) 


with A9 in radians given in Eq. ®. 

In Eqs. (ITTll . lIT^ and Ca, kE, kL, ^A, kR, a and f3 
are parameters to be determined by fitting the results of 
Monte Carlo simulations of radiopulse emission in differ¬ 
ent media to the parameterizations. For this purpose, 
we have performed fits to the frequency spectrum at the 
Cherenkov angle and away from it, as well as to the angu¬ 
lar distribution of the field around the Cherenkov angle, 
in ice, salt and the regolith using Eqs. ((TT|) . inj) and 
m, and the numerical values of p, Xq, Rm, Ec and n in 
Table E Table UTI lists the results of these fits for the dif¬ 
ferent media considered in this work. The fits are shown 
in Figs. 0 and 0 

As shown in Table QH it is remarkable that the nu¬ 
merical values of the normalization constants are weakly 
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TABLE II: Numerical values of ks (in V cm“^ MHz“^), fci,, 
kA, kn, a and P after fitting the results of Monte Carlo sim¬ 
ulations to Eqs. iiTTt . jnj and im in different media. 


Medium 

ks 

kL 

kA 

ka 

a 

P 

Ice 

4.79 10"^® 

23.80 

18.33 

1.42 

1.32 

2.25 

Salt 

3.20 10"^® 

21.12 

14.95 

1.33 

1.27 

2.60 

Moon 

3.30 10-1® 

23.74 

17.78 

1.41 

1.23 

2.70 


dependent on the medium. This is reflecting that most 
of the dependence on the medium is already accounted 
for in the unified parameterizations through the param¬ 
eters p, Xq, Rm, Ec and n. This is true at ^ 40% level 
in kE, and at ^ 15% or less in the other normalization 
constants, which is remarkably good given the simplic¬ 
ity of our toy model of shower development in which the 
parameterizations are based. 

The largest departure from being medium independent 
appears in The reason for this is that we made the 
simplification that particles in the shower follow straight 
lines parallel to the shower axis, while it is well known 
that they deviate from the shower axis by an average 
angle which is medium dependent. This dependence is 
not included in Eq. cni. As a result, the projection of the 
tracks in the perpendicular to the direction of observation 
is not simply accounted for by a factor sin0, but some 
dependence on the deviation angle should be included 
in Eq. lig. The small departures from being medium 
independent in the other normalization constants can be 
attributed to the fact that an electromagnetic shower is 
better represented by a “cone” with medium-dependent 
opening angle, rather than by a “cylinder” as depicted 
in Fig. 121 This introduces an extra dependence on the 
medium which is not accounted for in Eqs.0 and O- 

Finally, we are assuming in Eq. m that the electric 
held scales with the tracklength produced by particles 
with energy above the critical energy. However the elec¬ 
tric field is known to scale with the excess charge track- 
length, which is mainly due to Compton scattering dom¬ 
inating at energies below the critical energy. Using the 
total tracklength is however not a bad approximation be¬ 
cause it scales in roughly the same way with medium pa¬ 
rameters as the excess tracklength. This is confirmed by 
our GEANT4 simulations which predict a ratio of total 
to excess tracklengths with almost no dependence on the 
medium, of the order of 3.07, 3.05 and 2.98 in ice, salt 
and the lunar regolith, respectively. 

C. Discussion: Range of applicability of the 
parameterizations 

Eq. (HU together with the numerical values of the nor¬ 
malization constants kE, fc/j and a, describes accurately 
the electric field spectrum at the Cherenkov angle for 
frequencies below 10 GHz. This is well above the fre¬ 


quencies of interest to experiments using ice, salt and the 
lunar regolith as targets for neutrino detection. The pa¬ 
rameterization is also approximately valid at all energies, 
because the lateral dimension of the shower ~ responsi¬ 
ble for the turnover frequency at the Gherenkov angle - 
does not change much with shower energy as Monte Garlo 
simulations of high energy showers in ice have confirmed 

mm. 

Away from the Cherenkov angle Eq. together 

with the numerical values of the normalization constants 
kE, kE and /3, describe accurately the electric field spec¬ 
trum up to about twice the frequency at which the 
turnover in the spectrum occurs. However, contrary to 
Eq. ([TTll . Eq. llT^ is expected to fail at energies above 
the energy (Et.pm) at which the Landau-Pomeranchuk- 
Migdal effect starts to affect the longitudinal shower 
development of the shower. The LPM effect is known to 
stretch the longitudinal development of an electromag¬ 
netic shower, and to change its effective longitudinal di¬ 
mension. The LPM effect also produces a decrease of 
the turnover frequency in the spectrum as can be easily 
understood from Eq. o and was shown in [m- In the 
last column of Table HI we list Eepm in ice, salt and the 
lunar regolith, using the definition of Eepm in [4^. 

The angular distribution of the electric field around the 
Gherenkov angle in Eqs. m and m, is also determined 
by the longitudinal development of the shower, and as a 
consequence it is also affected by the LPM effect. The 
parameterization in Eqs. and m is only valid below 
It describes - with an accuracy of ~ 10% - the 
angular distribution of the electric held for frequencies 
above ^ 100 MHz, and in the angular range where the 
pulse falls by a factor ^ 2-3 with respect to its value at 
the peak. 

Our parameterizations are also expected to be a good 
approximation to the coherent Cherenkov radio emis¬ 
sion in hadronic showers even at energies above Elpm- 
Hadronic showers are known to be less affected by the 
LPM effect m , and their efficiency of emission of coher¬ 
ent Cherenkov radiation at high energy has been shown 
to be similar to that of electromagnetic showers ji^l • H®" 
sides, at high energy, a large fraction of the energy of 
the hadronic shower goes into its electromagnetic com¬ 
ponent |37| . Our parameterizations are then expected to 
be applicable in estimates of radio emission in most ex¬ 
periments looking for extremely high energy neutrinos. 
These experiments are currently concentrating most of 
their efforts on the detection of the hadronic showers 
produced in charged current interactions, or in the 
neutral current interaction of any neutrino flavor. The 
observation of the electromagnetic shower produced in 
a Ve charged current interaction is expected to be very 
difficult because the LPM effect shrinks the angular dis¬ 
tribution of the electric field reducing dramatically the 
available solid angle for detection [i^. 

The important issues in the last two paragraphs above, 
call for simulations of electromagnetic and hadronic 
showers with energies above Elpm in different dense me- 
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dia. We defer this work to a future paper [i^ . 
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